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Abstract — We consider the problem of energy-efficient trans- 
mission in cooperative multihop wireless networks. Although the 
performance gains of cooperative approaches are well known, 
the combinatorial nature of these schemes makes it difficult to 
design efficient polynomial-time algorithms for deciding which 
nodes should take part in cooperation, and when and with 
what power they should transmit. In this work, we tackle 
this problem in memoryless networks with or without delay 
constraints, i.e., quality of service guarantee. We analyze a wide 
class of setups, including unicast, multi-cast, and broadcast, and 
two main cooperative approaches, namely: energy accumulation 
(EA) and mutual information accumulation (MIA). We provide 
a generalized algorithmic formulation of the problem that en- 
compasses all those cases. We investigate the similarities and 
differences of EA and MIA in our generalized formulation. We 
prove that the broadcast and multicast problems are, in general, 
not only NP hard but also o(log(n)) inapproximable. We break 
these problems into three parts: ordering, scheduling and power 
control, and propose a novel algorithm that, given an ordering, 
can optimally solve the joint power allocation and scheduling 
problems simultaneously in polynomial time. We further show 
empirically that this algorithm used in conjunction with an 
ordering derived heuristically using the Dijkstra's shortest path 
algorithm yields near-optimal performance in typical settings. For 
the unicast case, we prove that although the problem remains NP 
complete with MIA, it can be solved optimally and in polynomial 
time when EA is used. We further use our algorithm to study 
numerically the trade-off between delay and power-efficiency in 
cooperative broadcast and compare the performance of EA vs 
MIA as well as the performance of our cooperative algorithm 
with a smart non-cooperative algorithm in a broadcast setting. 



I. INTRODUCTION 

In a wireless network, a transmit signal intended for one 
node is received not only by that node but also by other 
nodes. In a traditional point-to-point system, where there 
is only one intended recipient, this innate property of the 
wireless propagation channel can be a drawback, as the signal 
constitutes undesired interference in all nodes but the intended 
recipient. However, this effect also implies that a packet 
can be transmitted to multiple nodes simultaneously without 
additional energy expenditure. Exploiting this broadcast ad- 
vantage, broadcast, multicast and multihop unicast systems 
can be designed to work cooperatively and thereby achieve 
potential performance gains. As such, cooperative transmission 
in wireless networks has attracted a lot of interest not only 
from the research community in recent years [1], [2], [4], [5], 
[6], [7], [9] but also from industry in the form of first practical 
cooperative mobile ad-hoc network systems [14]. 
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We focus on the problem of cooperative transmission in 
this work, where a single node is sending a packet to ei- 
ther the entire network (broadcast), a single destination node 
(unicast) or more than one destination node (multicast), in 
a multihop wireless network. Other nodes in the network, 
that are neither the source nor the destination, may act as 
relays to help pass on the message through multiple hops. 
The transmission is completed when all the destination nodes 
have successfully received the message. We particularly focus 
on the case where there is a delay constraint, whereby the 
destination node(s) should receive the message within the 
delay constraint, however, we also discuss how our results 
apply to the unconstrained case. 

A key problem in such cooperative networks is routing 
and resource allocation, i.e., the question which nodes should 
participate in the transmission of data, and when, and with 
how much power, they should be transmitting. The situation 
is further complicated by the fact that the routing and resource 
allocation depends on the type of cooperation and other 
details of the transmission/reception strategies of the nodes. 
We consider in this paper a time-slotted system in which the 
nodes that have received and decoded the packet are allowed 
to re-transmit it in future slots. During reception, nodes add 
up the signal power (energy accumulation, EA) or the mutual 
information (mutual information accumulation, MIA) received 
from multiple sources. EA, which has been discussed in 
prior work [2], [5], [4], [6], can be implemented by using 
maximal ratio combining (MRC) of orthogonal signals from 
source nodes that use orthogonal time/frequency channels, or 
spreading codes, or distributed space-time codes. MIA can 
be achieved using rateless codes [9], [10]. Although these 
techniques are often treated separately in the literature, we 
shall see how our formulation of the problem encompasses 
both approaches and allows many of the results to be extended 
to both. 

We furthermore assume that the nodes are memoryless, i.e., 
accumulation at the receiver is restricted to transmissions from 
multiple nodes in the present time slot, while signals from 
previous timeslots are discarded. This assumption is justified 
by the limited storage capability of nodes in ad-hoc networks, 
as well as the additional energy consumption nodes have to 
expand in order to stay in an active reception mode when they 
overhear weak signals in preceding timeslots. Note that much 
of the literature cited above has used the assumption of nodes 
with memory, so that their results are not directly comparable 
to ours. 

A key tradeoff is between the total energy consumption 1 

'As we consider fixed time slot durations, we use the words energy and 
power interchangeably in this paper 
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and the total delay measured in terms of the number of slots 
needed for all destination nodes in the network to receive the 
message. At one extreme, if we wish to minimize delay, each 
transmitting node should transmit at the highest power possible 
so that the maximum number of receivers can decode the mes- 
sage at each step (indeed, if there is no power constraint, then 
the source node could transmit at a sufficiently high power to 
reach all destination nodes in the first slot itself). On the other 
hand, reducing transmit power levels to save energy, may result 
in fewer nodes decoding the signal at each step, and therefore 
in a longer time to complete the transmission. We therefore 
formulate the problem of performing this transmission in such 
a way that the total transmission energy over all transmitting 
nodes is minimized, while meeting a desired delay constraint 
on the maximum number of slots that may be used to complete 
the transmission. The design variable in this problem is to 
decide which nodes should transmit, when, and with what 
power. 

The key contributions of our work are as follows: 

> We formulate the problem of minimum energy transmis- 
sion in cooperative networks. Although the prior literature 
have focused on either EA ([2], [6]) or MIA ([9], [10]) 
and have treated them separately, our generalized formu- 
lation can treat both methods as variations of the same 
problem. 

• Our formulation of delay-constrained minimum energy 
broadcast in cooperative networks, goes beyond the prior 
work in the literature on cooperative broadcast which 
has focused either on minimizing energy without delay 
constraints [2], [6], or on delay analysis without energy 
minimization [8]. Our extended problem formulation 
allows us to expose and investigate the energy-delay 
tradeoffs inherent in cooperative networking. 

• We not only prove that the delay constrained minimum 
energy cooperative broadcast (DMECB) and multicast 
(DMECM) problems are NP-complete, but also that they 
are o(log(n))-inapproximable (i.e., unless P = NP, it is 
not possible to develop a polynomial time algorithm for 
this problem that can obtain a solution that is strictly 
better than a logarithmic-factor of the optimum in all 
cases). We are not aware of prior work on cooperative 
broadcast or multicast that shows such inapproximability 
results. 

• We show that the delay constrained minimum energy 
cooperative unicast (DMECU) problem is solvable in 
polynomial time using EA but is NP-complete using 
MIA. We are unaware of any hardness results on unicast 
approaches using mutual information accumulation. 

• For the cases where we prove the transmission problem to 
be hard, we are able to show that for any given ordering 
of the transmissions (which dictates that a node later in 
the ordering may not transmit before the nodes earlier in 
the ordering have decoded successfully), then the problem 
of joint scheduling and power allocation can in fact be 
solved optimally in polynomial time using a combination 
of Dynamic Programming for the scheduling and Convex 
Programming for the power allocation. 



• For small network instances, we compute the optimal 
solution through exhaustive search, and show empirically 
through simulations that our proposed joint scheduling 
and power control method works near-optimally in typical 
cases when used in conjunction with an ordering provided 
by the Dijkstra tree construction. 

• We also show through simulations the delay-energy 
tradeoffs and minimum energy performance for larger 
networks and demonstrate the significant improvements 
that can be achieved by our solution compared to non- 
cooperative broadcast. We further compare the perfor- 
mance of our proposed broadcast algorithm under MIA 
and EA approaches. 

• For DMECB where EA is used, we present a reduction 
that would allow for a polynomial time algorithm for 
the joint ordering-scheduling-power control problem that 
is provably guaranteed to offer a 0(n e ) approximation, 
for any e > 0. This algorithm is based on the current 
best-known algorithm for the bounded diameter directed 
Steiner tree problem [19]. Using the same reduction, we 
can also get an approximation factor of 0(T log 2 (n)) 
for a fixed delay constraint T. Given that DMECB is 
o(log(n)) inapproximable for any T > 2, this provides a 
fairly tight approximation, especially for small T. 

This paper makes several contributions that significantly 
enhance our understanding of complexity and algorithm de- 
sign for cooperative transmission in wireless networks, in 
the context of energy-delay tradeoffs. The summary of the 
algorithmic results developed in this paper are presented in 
Tables I and II. It is worth noticing that although in this work 
we have considered a memoryless setting, all the negative 
results presented in Table I extend to the case where there 
is no memory. 

The rest of the paper is organized as follows: section II 
places our work in the context of prior related work. Section III 
describes the system model. The problem formulation is 
presented in section IV. In section V, we establish hardness 
results. The main positive result for the hard problems de- 
scribed in V is presented in section VI, where a polynomial 
time algorithm for optimum delay constrained scheduling and 
power allocation is presented, for the case when the ordering 
is given. In section VII, we provide a polynomial time optimal 
algorithm for the unicast case where EA is used. Section VIII 
provides approximation results for DMECB when EA is used. 
Simulation results are presented in section IX, where we 
suggest and evaluate several heuristics for the ordering. We 
present concluding comments and directions for future work 
in section X, along with two tables summarizing the results. 

II. RELATED WORK 

EA and MIA are two of the main approaches 2 for cooper- 
ative communications in wireless network . 

In EA [2], [5], [4], [6], a receiver can recover the original 
packet so long as the total received energy from multiple 
sources or successive transmissions exceeds a given threshold. 

2 Other cooperation schemes include distributed beamforming [1] and coded 
cooperation [3], which will not be further considered in this paper. 
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NEGATIVE RESULTS 


Energy Accumulation 


Mutual Information Accumulation 


Delay Constraint (T) 


Unconstrained 


Delay Constraint (T) 


Unconstrained 


Broadcast 


o(log(ra)) inapproximable for T > 3 


NP-complete 


o(log(n)) inapproximable for T > 3 


NP-complete 


Multicast 


o(log(n)) inapproximable for T > 3 


NP-complete 


o(log(n)) inapproximable for T > 3 


NP-complete 


Unicast 


Polynomial time 


Polynomial time 


NP-complete for T > 4 





TABLE I: Summary of the algorithmic negative results. 



POSITIVE RESULTS 


Energy Accumulation 


Mutual Information Accumulation 


Delay Constraint (T) 


Unconstrained 


Delay Constraint (T) 


Unconstrained 


Broadcast 


•0{n e ) for e > 
•0(T log 2 (n)), for fixed 
T 

• Polynomial time given 
ordering (DMCT_go) 

• Polynomial time algo- 
rithm for T = 1 , 2 


Polynomial time given 
ordering 

ID dynamic program 
given in (7) 


• Polynomial time given 
ordering (DMCT_go) 

• Polynomial time algo- 
rithm for T = 1, 2 


Polynomial time given 
ordering 

ID dynamic program 
given in (7) 


Multicast 


Polynomial time given or- 
dering (DMCT_go) 


Polynomial time given or- 
dering (DMCT_go) 


Polynomial time given or- 
dering (DMCT_go) 


Polynomial time given or- 
dering (DMCT_go) 


Unicast 


Polynomial time 


Polynomial time 


• Polynomial time given 
ordering (DMCT go) 

T > 4 

• Polynomial time 
algorithm for T = 1, 2 





TABLE II: Summary of the algorithmic positive results. 



Such an approach can be implemented using maximal ratio 
combining of orthogonal signals from multiple sources, e.g., 
through a Rake receiver in CDMA or distributed space-time 
codes. It has been shown that one can achieve significant sav- 
ing in energy and/or transmission time when using an energy 
accumulation protocol, compared to traditional protocols [2], 
[6], [8]. If energy accumulation is achieved by transmitting 
the exact same packet either from different relays or through 
successive re-transmissions, the scheme is shown to achieve 
capacity in an asymptotically wideband regime [2]. We note 
that recently a commercially developed cooperative mobile 
ad hoc network system has been developed which utilizes a 
pragmatic cooperation method requiring minimal information 
exchange, based on a combination of phase dithering and turbo 
codes [13], [14]. It is shown in [13] that the performance of 
this pragmatic scheme is close to that of an ideal energy- 
accumulation approach based on space-time coding. 

In MIA [9], [10], the receiving node accumulates mutual 
information for a packet from multiple transmissions until it 
can be decoded successfully. In practice, this can be achieved 
using rateless codes [9], [10]. The two schemes have been 
shown to be equivalent at low signal-to-noise ratios (SNRs), 
while in high SNR regimes the latter is shown to have superior 
performance [9], [10]. 

Many network protocols in mobile ad hoc and sensor 
networks need to operate in broadcast mode to disseminate 
certain control messages to the entire network (for instance, 
to initiate route requests, or to propagate a query). The subject 
of broadcast transmission in multi-hop wireless networks has 
attracted a lot of attention from the research community in both 
non-cooperative [15], [16], [18] and cooperative settings [2], 
[4], [6], [8], [7]. For traditional non-cooperative wireless 
networks, Cagalj et al. [18] showed that the problem of 
minimum energy broadcast is NP-hard. In [4], Mergen et 
al. show through a continuum analysis the existence of a 
phase transition in the behavior of cooperative broadcast: if the 



decoding threshold is below a critical value then the broadcast 
is successful, else only a fraction of the network is reached. 
In [6], Mergen and Scaglione, show that the problem of 
scheduling and power control for minimum energy broadcast 
is tractably solvable for highly dense (continuum) networks 
and show the gains obtained with respect to noncooperative 
broadcast. In [8], we examined the delay performance of 
cooperative broadcast and show that cooperation can result 
in extremely fast message propagation, scaling logarithmically 
with respect to the network diameter, unlike the linear scaling 
for non-cooperative broadcast. 

The work by Marie and Yates [2] is closest in spirit to 
our work. They too address the minimum-energy cooperative 
broadcast problem where nodes use EA. However, in their 
work the system has memory, in that the nodes can save soft 
information from all previous transmissions throughout time 
and use it to decode data later on. They prove that the problem 
is NP-complete in this case. In their setting, because of the 
memory, it suffices to have each transmitter transmit only once; 
therefore there is no distinction between ordering and schedul- 
ing. This is no longer true in our memoryless setting where 
the energy from past transmissions cannot be accumulated. A 
further key distinction in our work is that we consider delay 
constraints, whereas [2] focuses only on the minimum energy 
cooperative solution without delay constraints; additionally we 
consider MIA, which is not addressed in [2]. 

One prior work that discusses the power-delay tradeoff in a 
cooperative setting is [17]; however, the focus of that work is 
on space-time codes used for unicast, not broadcast. In earlier 
work [12], we had considered this tradeoff in a broadcast 
cooperative setting using EA and conjectured that many of 
the results would extend to MIA but the investigation of that 
conjecture had remained an open problem. In this paper, we are 
addressing that open problem, as well as considering unicast 
and multicast settings - which were not addressed in [12]. 

MIA is most notably discussed in [9], [10]. However, the 
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discussions are focused on unicast routing. The algorithms 
presented are heuristics, the performance of which are verified 
via simulations. We are not aware of any hardness results on 
cooperative MIA or any discussions on the use of MIA in a 
general cooperative broadcast setting. 

III. SYSTEM MODEL 

We consider a wireless network with n nodes. Radio propa- 
gation is modeled by a given symmetric n by n static channel 
matrix, H = {hij}, representing the (power) gain on the 
channel between each pair of nodes i and j. Time is assumed 
to be discretized into fixed-duration slots; without loss of gen- 
erality we assume unit slot durations. We assume cooperative 
communication in the receivers, encompassing two scenarios: 
EA and MIA. Only a single message is transmitted through 
the network. 

In EA, the received power at a given receiver in a specific 
timeslot is sum of the powers received from the transmitters 
that are active during that slot. As described in [4], [6], [5], 
this kind of additive received power can be achieved via 
maximal ratio combining under different scenarios including 
transmission using TDMA, FDMA channels, as well as with 
CDMA spreading codes and space-time codes. 

MIA can be implemented using rateless codes and decoders 
at receivers, as described in [9], [10], [11]. With proper design 
(e.g., different spreading codes), information streams from 
different relay nodes can be distinguished, and the mutual 
information of signals transmitted by different nodes can be 
accumulated. We consider a per-node bandwidth constraint and 
dynamic power allocation. 

We assume appropriate coding is used so that each receiving 
node can decode the message so long as its accumulated 
received mutual information exceeds a given threshold 9 that 
represents the bandwidth-normalized entropy of the informa- 
tion codeword in nats/Hz. Furthermore, all nodes are assumed 
to operate in half-duplex mode, i.e. they cannot transmit and 
receive simultaneously. If used in transmission, the nodes 
operate based on a decode and forward protocol. Therefore, 
they are not allowed to take part in transmission until they 
have fully decoded their message. 

Assuming the noise power is the same at all receivers, we 
can assume without loss of generality the noise power to be 
normalized to unity so that the transmit power attenuated 
by the channel becomes equivalent to the signal to noise 
ratio (SNR). As mentioned in Sec. I, we assume a memory- 
less model in which nodes do not accumulate energy or 
information from transmissions occurred in previous time 
slots. 

IV. PROBLEM FORMULATION 

In this section we provide a generalized formulation for the 
delay constrained minimum-energy cooperative transmission 
(DMECT) problem in the setting described in Sec. III. 

We assume that the transmission begins from a single source 
node. The aim is to get the message to all the nodes in a 
destination set V, with the minimum possible total energy, 
within a time T (which can take on any value from 1 to 



n — 1). Every node in the network is allowed to cooperate 
in the transmission, so long as they have already decoded the 
message. The problem now becomes: which nodes should take 
part in the cooperation, when and with what power should they 
transmit to achieve this aim while meeting the constraints and 
incurring minimum total transmission power. 

Recalling the memoryless assumption, the condition for 
successful decoding at some receiver node r at time t when a 
set of nodes S(t) is transmitting packets, with transmit power 
p st ,Vs e S(t) is: 

Vrt > (1) 

with y rt being the mutual information accumulated by node r 
at time t. Let xa be an indicator binary variable that indicates 
whether or not node i is allowed to transmit at time t. In other 
words, we define xu to be 1, if node i is allowed to transmit 
at time t (i.e. has decoded the message by the beginning of 
time slot t as per equation (1)), and otherwise. Let pu be 
the transmit power for each node i at each time t. Without 
loss of generality, the source node is assigned node index 1. 

The DMECT problem can then be formalized as a combi- 
natorial optimization problem: 

min P total = Ylt=i S"=i Pit ( 2 ) 

s.t. 1. p it > 0, Vi,Vi 

2. x iT+1 > 1, Mi e V 

3. Xit+i < \yu +x it , Vi,Vt 

4. x lt = 1, Vi 

5. x a = 0, Vi ^ 1 

6. x lt e{0,l} 

where, for the energy accumulation (EA) case 3 : 

y it = log 1 + ^2 Pstx s th sl (3) 

and for mutual information accumulation (MIA) case: 

Vit = ^2 \og(l+p st x st h sl ) (4) 

ses(t) 

Constraint 2 ensures that every nodes in the destination set 
successfully decodes the message within the time constraint 
T, constraint 3 ensures that a node cannot transmit unless it 
has already received the message while simultaneously making 
sure that a node that has decoded the message in previous time 
slots will not be prevented from transmitting in future time 
slots (if it wants to transmit), constraint 4 assigns the source 
node, and all other constraints are self-explanatory. 

In general, there are three variations of this problem, based 
on the size of the destination set: 

• Delay constrained minimum energy cooperative unicast 
(DMECU): where the set V includes a single destination 
node. 

3 Notice that because of the monotonicity of the log function, y it > 9 in 
this case is equivalent to X)seS(t) Pstx s th s i > e e — 1 
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• Delay constrained minimum energy cooperative multicast 
(DMECM): where the set T> includes more than one 
destination node. 

• Delay constrained minimum energy cooperative broad- 
cast (DMECB): where the set V includes all the nodes 
apart from the source node. 

The decision version of these problem, can be defined 
correspondingly as follows: "Given some power bound C, 
does there exist an allocation of powers, p it , satisfying the 
constraints in (2) such that P to tai < C?" An instance of this 
decision problem is defined by giving the symmetric n x n 
matrix H, with a designated source node (vertex), a destination 
set V, a delay bound T, and a power bound C. 

Notice that assigning T > n, in the above formulation, 
results in the problem definition in the case where there is 
no delay constraint. Note also that a requirement for per- 
node maximum power can be trivially added to the above 
formulation as additional constraint; we have left that out for 
simplicity. Should the maximum power be added, it should be 
large enough to ensure a feasible solution exists for the given 
connectivity and delay constraint. 

V. HARDNESS 
In this section we prove that finding an optimal solution 
for DMECB and DMECM problems is not only NP-hard but 
also o(log(n)) inapproximable i.e., finding any polynomial 
time algorithm that approximates the optimal solutions within 
a factor of o(log(n)) is also NP-hard. We show this by 
demonstrating that any instance of the set cover problem can 
be reduced to an instance of DMECB (and by extension 
DMECM). We further prove NP-completeness for DMECU 
when MIA is used; note that DMECU with EA will be treated 
in Sec. VII. 

A. Set Cover Problem 

The set cover problem is a classical problem in computer 
science [20]. It is stated as follows: Given a universe U of n 
elements and a collection of subsets of U, S = S\, S2, ---Sk, 
find a minimum subcollection of S that covers all elements of 
U. This problem is NP-complete and was shown, in [21], to 
be o(log(n)) inapproximable. 

The set cover problem can be thought of as a bipartite graph 
G(V,E), with \V\ = k + n, representing the k sets and n 
elements and the edges are used to connect each set to its 
elements. This is shown in Figure 1 (a), where we assign a 
vertex for each set in the top part of the graph, and assign a 
vertex for each element in the bottom part of the graph. We 
connect each set to its elements using an edge. One can think 
of each vertex in this graph as a node in a network, in which 
edges exist between any pair of nodes for which hij > 0, 
and the edges are labeled with a weight Wij that corresponds 
to the transmit power needed at node i to exceed a threshold 
of 9 at the receiver j, in a single time slot if i was the only 
transmitter. Given an instance, G, of the set cover problem, 
the optimal solution to the set cover problem, OPT sc , would 
find the minimum subset of vertices in the top part of the 
graph, so that their transmission of a message can broadcast 
the message to all the vertices in the bottom part of the graph. 



B. Inapproximablity of DMECB 

Given an instance, G, of the set cover problem, with k sets 
and n elements, let us construct a new graph G' as follows: 
Assign a root node r, which is the source with the message at 
the starting time, call this level 0. Include k nodes in level 1, 
representing the k sets in the set cover problem, all connected 
to the root node, as shown in Figure 1 (b). This is followed by 
the bipartite graph of G, which makes up level 2 and 3 of G'. 
Connect each of the k nodes in level 2 to their representative 
in level 1 and to all the other nodes in level 2. Notice the nodes 
in level 2 are also connected to their elements in level 3 of 
the graph, as shown in the Figure. We make all the weight on 
the edges arbitrarily small (say 1), with the exception of the 
edges in between the nodes in level 1 and 2. We make those 
edges to be sufficiently large, say M, to be specified later. 



r 



G 


G' 
M 


-A 

M 




D 

M 


Level 
Level 1 














Level 2 
Level 3 







Fig. 1: Construction of G' for a given G, notice that not all 
the edges are shown (for clarity). 



Assume the the weight on the edges represent the power 
needed for the message to be transmitted across that edge. If 
we were to run the optimal DMECB algorithm on G' with 
T = 3 the algorithm would have to act as follows, to be able 
to cover all the nodes in the given time frame: 
Step 1: Root transmits with power 1, turning on all its k 
neighbors on level 1. 

Step 2: The algorithm picks a subset of the k nodes on level 1 
to transmit the message. This subset must be chosen to be as 
small as possible, given the large weight they have to endure to 
pass on the message on to the bipartite graph, and the fact that 
DMECB is trying to minimize the total weight. Yet it has to 
be large enough so that when the nodes in level 2 transmit, all 
the nodes in level 3 would receive the message. The optimal 
algorithm must be able to find such a subset. 
Step 3: The nodes that receive the message in level 2 transmit 
the message in this step, turning on all the nodes in level 3 of 
the graph, as well as all the nodes in level 2 of the graph that 
were not selected for transmission, thus covering the whole 
graph. 

Let us call the solution 4 of this optimal algorithm 
OPT DMEC b- Then the following two lemmas with respect 
to the above construction of G and G' hold: 

Lemma V.l OPT DM ecb < M.OPT sc + 1 + OPT sc 

4 Minimum energy needed for transmission. 
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Proof: Consider an instance of SC (with graph G), 
whose optimal solution is OPTsc- Construct a graph G', as 
explained and run the DMECB algorithm to get OPT D mecb- 
The above inequality holds by construction of the graph. ■ 

Lemma V.2 OPT sc < OPT ™ E ° B 

Proof: Consider an instance of DMECB on G' and its 
optimal solution OPTdmecb for delay T = 3. Notice that 
if T > 3, we add additional single nodes (as virtual roots) 
to reduce the problem to the case where T = 3. Looking 
at G', we observe that to meet the delay constraint, by end 
of step i, at least one node in level i must have heard the 
message - else it is impossible to get the message through 
to the rest of the levels in the time frame left. Let's say the 
root is on level 0. Consider the subset of level 1 that has 
come on at the end of time 1, Si, and from level 2 consider 
the set, s 2 , that came on at the end of time step 2. We now 
want to show that s 2 is a feasible solution for set cover. To 
do so, we make the following two claims: Claim 1: Nodes 
responsible for turning on s 2 must be a subset of Si . Claim 2: 
S2 is a feasible solution to set cover. Claim 1 holds because 
only nodes that have received the message by the end of time 
1 can transmit the message at time 2. Not all of them might 
transmit though, so s 2 is a subset of corresponding nodes in 
Si. Claim 2 is true because if there exists an element in level 
3 that is not a corresponding node to any node in s 2 , it cannot 
decode by T = 3. Therefore, s 2 , is a feasible solution to set 
cover. OPTdmecb must spend at least M for each element 
of s 2 to come on, so OPT sc < Q^Mmm _ a 

Theorem V.l The DMECB problem is o(log(n)) inapprox- 
imable, for T > 3. 

Proof: For an instance of the set cover problem, with k 
being the total number of sets, lemma V.l can be re-written 
as OPT D mecb < M.OPTsc + l + k. We also know by 
lemma V.2 that OPT sc < OPTd ^ ecb . Therefore, for a suffi- 



ciently large M, we can write OPTsc 



OPT, 



DMECB 



M 



■0(1). 



Therefore, the reduction used in construction of the graph G' 
preserves the approximation factor. That is, if one can find an 
a-approximation for DMECB, by extension there must exist 
an a-approximation for set cover. We know, by [21], that the 
set cover problem is o(log(n)) inapproximable, thus DMECB 
must be o(log(n)) inapproximable. In other words, finding a 
polynomial time approximation algorithm that approximates 
OPT D mecb with a factor of o(log(ra)) is NP-hard. ■ 
The DMECB problem can be solved in polynomial time 
for cases when T < 3. The optimal algorithm for T = 1 
is trivial and an optimal polynomial algorithm for T = 2 
is discussed in section VI 5 . It is also trivial to verify the 
feasibility of a given power allocation, and verify whether 
or not it satisfies the decision version of DMECB given in 
section IV. Therefore, the problem belongs to the class of NP. 
Notice that the inapproximability result, given by Theorem 

5 DMECT_go algorithm, discussed in section VI, along with an ordering 
based on channel gains from the source, provides an optimal polynomial time 
algorithm for DMECB for the case when T = 2. 



V.l, is stronger than, and implies, the NP-completeness result. 
It is also worth noticing that without any delay constraint (i.e. 
when T > n), the problem is still NP-complete and the proof 
can be obtained, using directed Hamiltonean path, following 
the approach in [2]. 

C. Inapproximablity of DMECM 

The proof of the following theorem, follows from Theorem 
V.l by noticing that broadcast can be thought of as a special 
case of multicast. 

Theorem V.2 The DMECM problem is o(log(n)) inapprox- 
imable, for T > 3. 

D. Hardness Results for DMECU 

In the unicast case, the hardness of the problem depends on 
whether we are using EA or MIA. In the former case, DMECU 
can be shown to be polynomially solvable and the algorithm 
for that is provided in section VII. In the remainder of this 
section, we discuss DMECU with MIA. 

Given an instance, G, of the set cover problem, with k sets 
and n elements, similar to that in section V-B, let us construct 
a new graph G' as follows: Assign a root node r, which is 
the source with the message at the starting time, call this 
level 0. Include k nodes in level 1, representing the k sets 
in the set cover problem, all connected to the root node with 
a small weight (say weight 1), as shown in Figure 2. This is 
followed by the bipartite graph of G, which makes up level 2 
and 3 of G'. Connect each of the k nodes in level 2 to their 
representative in level 1 with edge weights, of say W. Notice 
the nodes in level 2 are also connected to their elements in 
level 3 of the graph, as shown in the Figure, with low-weight 
edges. Add a single destination node d, in level 4 and connect 
all the nodes in level 3 to d. Let the channel between all nodes 
on level 3 and destination d be equal and of gain h. Therefore, 
the edge weight on the edges connecting the level 3 nodes to 
d, can be assigned to be M, where M is defined so that the 
following equality holds: log(l + Mh) = 9. 




Level Level 1 Level 2 Level 3 Level 4 



Fig. 2: Construction of G' for a given G in DMECU, notice 
that not all the edges are shown (for clarity). 

Assume that the weight on the edges represent the power 
needed for the message to be transmitted across that edge. If 
we were to run the optimal DMECU algorithm on G' with 
T = 4 the algorithm would have to act as follows, to be able 
to turn node d on within the given time frame: 
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Step 1: Root transmits with power 1, turning on all its k 
neighbors on level 1. 

Step 2: The algorithm picks a subset the k nodes on level 1 
to transmit the message to the nodes in level 2. 
Step 3: A subset of nodes that have received the message in 
level 2, transmit the message in this step, turning on a subset 
of nodes in level 3 of the graph. 

Step 4: A subset of nodes that have received the message in 
level 3, transmit the message in this step with sufficient power 
to turn on d. 

Let us call the solution of this optimal algorithm 

OPTdmecu ■ 

Theorem V.3 The DMECU problem, with MIA, is NP- 
complete for T > 4. 

Proof: Given an instance of G, we construct G' as above. 
Let us run DMECU on G' and call the optimal solution 
OPTdmecu for delay T = 4. Notice that if T > 4, we 
add additional single nodes (as virtual roots) to reduce the 
problem to the case where T = 4. Define p to satisfy the 
following: nlog(l + ph) = log(l + Mh), meaning p is the 
power required for nodes on level 3 to turn on d, if all of 
them were transmitting at the same time. Claim: OPTdmecu 
needs to use all the nodes in level 3 for transmission. This 
claim holds by contradiction, as follows: If all the nodes on 
level 3 are used for transmission, each node on that level 
must transmit with power p. Let's assume one of the nodes in 
that level is not used for transmission. Then the remaining 
nodes in level 3 need to transmit with power p', where 
nlog(l+ph) = (n— 1) log(l+p'h) = 9. Therefore, the ratio of 
the sum power needed with one fewer node transmitting to the 
case where all nodes in level 3 are transmitting can be written 

as {J ^ ]L = iI ^SfSr 1 ' for sufficientl y lar s e d > this 

ratio can become arbitrarily large. Therefore, for sufficiently 
large 9, the claim holds. Given the claim holds, we know that 
by definition OPTsc provides the optimal way (minimum 
energy) to turn on all the nodes in level 3 within the required 
time frame, therefore, for non-zero edge weights OPTdm ecu 
needs to optimally solve the set cover problem in step 2. ■ 
It is worth noticing that all the hardness results presented 
in this section extend to the case where there is no memory. 

VI. OPTIMAL Transmission GIVEN ORDERING 

In Section V, we proved NP-hardness for several variations 
of the DMECT problem (namely, the DMECB, DMECM and 
DMECU (with MIA), with the former two being o(log(n)) 
innaproximable). In this section, we break this NP-hard prob- 
lem into three subproblems, namely ordering, scheduling and 
power allocation, and we propose an optimal polynomial time 
algorithm for joint scheduling and power allocation when the 
ordering is given. We evaluate a heuristic for the ordering in 
Section IX. 

Definition An ordering, for a vector of n nodes, is an array 
of indices from 1 to n; any node that has decoded the message 
will only be allowed to retransmit when all nodes with smaller 



index have also decoded the message (and are thus allowed to 
take part in transmission). 

Given an ordering, what remains to be determined is which 
nodes should take part in transmission, how much power 
they should transmit with and at what time slots, such that 
minimum energy is consumed while delay constraints are 
satisfied. 

A. Instantaneous optimal power allocation 

If we know which nodes are transmitting the message and 
which nodes are receiving it, at any single time-slot, we can 
use a convex program (CP) to determine the optimal power 
allocation for that time slot. Consider an ordered vector of n 
nodes (1, k, i, n). Let us assume that by time slot t, 
node 1 to i have decoded the message and nodes i + 1 to n 
are to decode it during that time slot. At time instance t, the 
optimal instantaneous power allocation for a set of transmitting 
nodes (say S(t) = (k,...,i)) to turn on a set of receiving 
nodes nodes (say R(t) = (i + 1, n)) can be calculated by 
the following CP: 

min T, seS (t)Pst (5) 

S.t. p st > 0, Vs 

Vn >9,Vre R(t) 

We use the notation CP([{k...i},{i + l...n}],0,H) 
to refer to solution of the above CP. As a notation, 
CP([{x...y},{z...a}],9,H) = 0, if z > a. Notice that in 
the case where EA is used, this CP simply reduces to a linear 
program, using the manipulation highlighted in footnote 2 in 
section IV. 

B. Joint Scheduling and power allocation 

Knowing the instantaneous optimal power allocation given 
the set of senders and receivers at each time slot, all that 
remains to be done is to determine these sets at each time 
slot, in order to minimize the overall power while meeting the 
delay constraint. 

Let C (j, t) be the minimum energy needed to cover up to 
node j in t steps or less. We can calculate this, using the 
following algorithm: 

C(j, t) = min [C(k, t - 1) + CP({l...k}, {k + 9, H)] 

k£(l,..,j) 

(6) 

where C(k, 1) = CP(1, {2...k}, 9, H), C(l, t) = Vi. 

Thus, in DMECB, the total minimum cost for covering n 
nodes by time T can be found by calculating C(n,T). In 
DMECM, and DMECU (MIA), the same approach could be 
used, except for node n being replaced by the highest order 
destination in the former and by the destination order in the 
latter. 

A pseudocode for the algorithm is presented below: 
The complexity of the optimal scheduling and power allo- 
cation can be obtained by inspection of the above algorithm: it 
invokes at most 0(n 2 T) calls to the CP solver, each of which 
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Algorithm 1 Delay constrained minimum energy cooperative 
transmission, given an ordering (DMECT_go) 

1: INPUT: an ordered array of nodes of size n (where node 
i is the ith node in the array), T (delay), d (destination), 
H (channel), 6 (threshold). 

2: OUTPUT: C (cost matrix) 

3: Begin: 

4: for i := 2 to n do 

5: C(i,l):=CP([l,{2...i}],0,H); 

6: end for 

7: for t := 1 to T do 

8: C(l,i):=0; 

9: end for 

10: for t := 2 to T do 

11: for i := 1 to do 

12: for A: := 1 to i do 

13: := C(k,t - 1) + CP([{l...fc},{fc + 

l...i}],9,H); 

14: end for 

15: C{i, t) :~ minx 

16: end for 

17: end for 



takes polynomial time. Hence the DMECT_go algorithm that 
does joint scheduling and power-control is a polynomial time 
algorithm. 

Note that the delay constraint T can be made sufficiently 
large (> n), or removed entirely from the formulation, to 
cover the case of no delay constraints. In that case the two- 
dimensional dynamic program proposed in (6), reduces to a 
one-dimensional dynamic program: 

C(n) = min [C(jfe) + CP({l...k}, {k + l...n}), 9, H] (7) 

k 

where C(n) is the minimum cost of covering node n using 
our cooperative memoryless approach, starting from node 1 
and C(l) = 0. 

VII. OPTIMAL UNICAST WITH ENERGY 
ACCUMULATION 

In this section we propose an optimal polynomial time 
algorithms for solving the unicast problem with EA. 

Theorem VII.l In DMECU with EA, there exists a solution 
consisting of a simple path between source and destination, 
which is optimum. 

Proof: Let us prove by induction: For delay T = 1, 
the claim is trivially true, as the optimal solution is a direct 
transmission from the source, s, to the given destination, d. 
For T > 1, we prove the claim by induction. Assume that 
the claim is true for T = k — 1. Pick any node in the 
network as the desired destination d. If the message can be 
transmitted from source s to d with minimum energy in a 
time frame less than k, then an optimal simple path exists by 
the induction assumption. So consider the case when it takes 
exactly T = k steps to turn on d. The system is memoryless, 



so d must decode by accumulating the energy transmitted 
from a set of nodes, v, at time k. This can be represented as 
^og(l+J2 v e\P v i k hdvi) > 0. We observe that there must exist 
a node v a G v whose channel to d is equal or better than all the 
other nodes in v. Therefore, given hd Vo > hd Vi , Vv^ G v — {v } 
then log(l + ^ ev 

Pvikhdv )>log(l+£_ )>6. 
In other words, if we add the power from all nodes in v and 
transmit instead from v a , our solution cannot be worse. v Q 
must have received the message by time k — 1, to be able to 
transmit the message to d at time k. We know by the induction 
assumption that the optimal simple path solution exists from 
source to any node to deliver the message within k — 1 time 
frame. Thus, for T = k, there exists a simple path solution 
between s and d, which is optimum. ■ 
Notice that the above theorem holds in the case where there 
is no delay constraint as well. The proof follows an straight- 
forward modification of the above proof and is omitted for 
brevity. 

Corollary VII.l The Dijkstra's shortest path algorithm pro- 
vides the optimal ordering in the case of minimum energy 
memoryless cooperative unicast, when there is no delay con- 
straint. 

Proof: We have already established that an optimal min- 
imum energy solution exists between source and destination, 
which is a simple path. The well-known Dijkstra's shortest 
path algorithm can find the minimum cost simple path be- 
tween source and destination. Therefore, Dijkstra's algorithm 
provides the optimal ordering. ■ 
Using theorem VII.l we know that the optimal unicast solution 
from source to any destination in DMECU (with EA) is given 
by a simple path. The cost paid by the optimal solution can 
be calculated using the following algorithm: Let C(i, t) be the 
minimum cost it takes for source node s to turn on i, possibly 
using relays, within at most t time slots. Then we can write: 

C{i,t)= min [C(k,t- 1) + w(k -> i)} (8) 

keNr{i) 

with C(s,t) = 0, for all t and C(i, 1) = w(s — > i), where 
Nr(i) is the set that contains i and its neighboring nodes that 
have a non-zero channel to i, w(k — > i) represents the power 
it takes for k to turn on i using direct transmission. Given that, 
the solution to OPTdmecu (with EA) is given by C(d,T). 
Computing this lower-bound incurs a running time of 0(n 3 ). 

The unicast case (with EA), with no delay constraint, is still 
polynomially solvable. Given Theorem VII.l and Corollary 
VII.l, the optimal solution is simply the weight of the shortest 
path given by the Dijkstra's algorithm. 

It is worth noticing that the crux of the difference between 
DMECU with EA and with MIA, that allows the former to be 
polynomially solvable, while the latter is NP-complete, lies in 
the optimality of single-node transmission. Namely, in the EA 
case, the multi-transmitters single-receiver case (multi-single) 
makes no sense as explained above and instead it is optimal to 
put the combined power into the best channel. This allows for 
the overall solution to be a simple path. However, in the MIA 
case, the many-to-one transmission case does in fact make 
sense. That is due to the property of the log function, creating 
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an effect similar to what we observe in water-filling, where it 
is best to transmit from the best channel up until some point, 
then from the second best channel and so forth. 

VIII. APPROXIMATION ALGORITHM FOR 
BROADCAST WITH ENERGY ACCUMULATION 

In Section V, we proved that DMECB is NP-complete and 
o(log(n)) inapproximable, therefore it is hard to approximate 
DMECB to a factor strictly better than log(n). It is of 
theoretical interest to know how close we can get to the 
optimal solution, using a polynomial-time algorithm. In this 
section we show that existing approximation algorithms for the 
bounded-diameter directed Steiner tree problem can be used 
to provide an 0(n e ) approximation for DMECB in the case 
where EA is used. We do so by proposing an approximation- 
preserving reduction to the directed Steiner tree problem. 

The Steiner tree problem is a classic problem in combinato- 
rial optimization [20]. We focus on a variation of this problem, 
namely bounded diameter directed Steiner tree, defined as 
follows. Given a directed weighted graph G(V, E), a specified 
root r G V, and a set of terminal nodes X C V {\X\ = n), 
the objective is to find the minimum cost arborescence rooted 
in r and spanning all vertices in X, subject to a maximum 
diameter T. Diameter refers to the maximum number of 
edges on any path in the tree. Notice that the tree may 
include vertices not in X as well, these are known as Steiner 
nodes. Directed Steiner tree problem is known to be NP- 
complete and 0(log(n)) inapproximable [20]. In [19], the 
authors give the first non-trivial approximation algorithms for 
Steiner tree problems and propose approximation algorithms 
that can achieve an approximation factor of 0(n e ) for any 
fixed e > in polynomial time. To the best of our knowledge 
this is currently the tightest approximation algorithm known 
for this problem. 

In order to reduce a given instance of the DMECB to an 
instance of the Steiner tree problem, we first restrict DMECB 
by not allowing many-transmitter- to-one-receiver (many- to- 
many) transmissions. Notice that in the proof of theorem VII. 1, 
we had established that many-to-one transmissions can be re- 
placed with one-to-one transmissions without loss of optimal- 
ity. Therefore, by not allowing many-to-many transmissions, 
we are left with one-to-one and one-to-many transmissions. 
We call this an integral version of DMECB, DMECB-int. The 
integrality gap of the weighted set cover problem is shown 
to be log(n) [20]; it is straightforward to extend that result to 
show that DMECB-int also loses a factor of log(n), compared 
to optimal DMECB. 

Consider an instance of DMECB-int, G(V, E), with (\V\ = 
n) and s G V being the source node. To reduce this problem to 
an instance of directed Steiner tree problem, let us construct a 
new graph G', consisting of n clusters, x', each corresponding 
to each node in G. Let each cluster be a bipartite graph, with 
n nodes on the left (marked as " — ") and n nodes on the 
right (marked as " + "), as shown in Figure 3. The " — " 
nodes are intra-connected within a cluster with edges of weight 
0. In each cluster, x' G G' corresponding to node v G G, 
the " + " and " — " nodes on each level, i, of the bipartite 



graph are connected to each other with an edge of weight 
Wi, representing the power needed by the corresponding node 
v G G to turn on its i closest neighbors. The i + node is then 
connected, with edges of weight 0, to all the " — " nodes in 
the corresponding neighbor clusters. We further add a single 
root node, r G G', and connect it via a zero- weight edge to 
all the " — " nodes in the cluster corresponding to s, x' s . We 
assign the root r and one desired " — " node from each cluster 
as terminal nodes and all other nodes in G' as Steiner nodes. 




Fig. 3: A simplified example of how clusters are constructed 
in G'. 

Let us look at an example of this construction, say node 
t>i G G, whose closest 3 neighbors are (v2,V4,ve)- We have 
an equivalent cluster x[ G G" corresponding to node v\. x[ 
has 2n nodes, arranged in n levels. The weight between the 
two nodes in say level 3 is equivalent to the power it takes for 
V\ to turn on (v 2 , V4, vq). Furthermore, the node 3+ in cluster 
x[ is connected to the " - " nodes in clusters (x' 2 , x' A , x' 6 ) with 
edges of weight 0. This construction allows us to find a way 
to allow vi to transmit with different power levels, without 
knowing what those powers might be in advance. We first add 
a single root node, r, and connect it via a zero-weight edge 
to all the " — " nodes in the cluster corresponding to s. 

Run the directed Steiner tree algorithm on G" to obtain a 
solution. The solution must choose at least one node from each 
cluster, to meet the mandatory terminal nodes requirement. 
Recall that each cluster in G' corresponds to a node in G and 
that multi-multi was not allowed. To convert the solution of 
the Steiner tree algorithm on G' to a solution of DMECB-int 
on G, we look at the parent of each cluster, which is a " + " 
node in another cluster. Let's say we want to see which node 
turns on vq by looking at G'. We look at the parent of x' 6 
and see that it's 3+ G x[. So in G, we figure out that v\ 
must transmit with enough power to turn on 3 of its closest 
neighbor (1V3), and it is as a result of this transmission that 
Vq comes on. Going through all the clusters and their parents, 
we can establish an ordering and transmission power for all 
the nodes that should take part in cooperation in G, and thus 
we have a solution for DMECB-int. 

Theorem VIII.l For DMECB problem with EA, an 0(n e ) 
approximation ratio can be achieved in polynomial time, for 
any fixed e > 0. 

Proof: As mentioned, the directed Steiner tree is 
o(log(n)) inapproximable, and the best approximation algo- 
rithm currently available [19] gives an 0(n e ) approximation 
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on the optimal solution. We had already lost 0(log(n)) to con- 
vert DMECB to its integral form. The approximation algorithm 
proposed in [19] can approximate the optimal integral solution 
within 0(n e ). Therefore, using the above reduction, and 
applying the directed Steiner tree approximation algorithm, 
we can approximate the optimal solution to DMECB within 
0(rf x log(ra)), which is equivalent to 0(n e ). ■ 
The running time of the Steiner tree approximation algo- 
rithm is a function of e, and the tighter the approximation, the 
worse the running time. Similarly, using the above mentioned 
reduction, the following result holds by directly applying the 
approximation algorithms in [19]. Detailed discussions of the 
algorithms in [19] are beyond the scope of this paper. 

Theorem VIII.2 For any fixed T > 0, there is an algorithm 
which runs in time and gives an 0(T log 2 (n)) approx- 

imation of the DMECB with EA. 

IX. PERFORMANCE EVALUATION 

For the simulations, we focus on the broadcast case. We 
consider a network of n nodes uniformly distributed on a 15 
by 15 square surface. The transmission starts from a node, ar- 
bitrarily located at the left center corner of the network (0,7). 
The channels between all nodes are static, with independent 
and exponentially distributed channel gains (corresponding to 
Rayleigh fading), where hij denotes the channel gain between 
node i and j. The mean value of the channel between two 
nodes, hij, is chosen to decay with the distance between 
the nodes, so that hij = d^ 1 , with dij being the distance 
between nodes i and j and r\ being the path loss exponent. 
The corresponding distribution for the channel gains is then 
given by 

f (fa) = =exp 

hij V hij J 

Based on the intuition developed in section VII, we use the 
Dijkstra's shortest path algorithm as our ordering heuristic. 
Simulations are repeated multiple times with the same node 
locations but different fading realizations and average values 
are shown in the graphs. Notice that the minimum power 
calculated by different algorithms, shown on the y-axes of the 
graphs in this section, are normalized by unit power (rendering 
it unit-less). The value of 6 is, arbitrarily, chosen to be log(2) 
throughout this section. 

In Figure 4, we calculate the optimal ordering by brute-force 
for a small number of nodes and compare the performance 
of our algorithm, which uses Dijkstra's shortest path-based 
ordering, with the optimal performance. The results for the 
broadcast case, with EA, is shown in this figure. As can 
be seen, Dijkstra's algorithm provides a good heuristic for 
ordering in this example and will be used throughout this 
section. Note that, as can be seen in the figure, although 
the problem was proved to be o(log(rc)) inapproximable, it is 
possible to achieve near-optimal results in polynomial time in 
certain practical settings where the network does not have any 
pathological properties. The inapproximablity results contain 
all possible (including pathological networks) scenarios. 




4 4.5 5 5.5 6 6.5 7 7.5 B 

Number of Nodes 



Fig. 4: Performance with optimal ordering vs Dijkstra's 
algorithm-based heuristic ordering. 

We next compare the performance of our cooperative al- 
gorithm with a smartly designed non-cooperative algorithm 
(both using EA). Notice that in our cooperative algorithm we 
make use of the wireless broadcast advantage (WBA), where 
transmission by one node can be received by multiple nodes 
and cooperative advantage, where a node can accumulate 
power from multiple transmitters. If an algorithm is using the 
WBA, but not the cooperative advantage, it can be thought 
of as an integral version of DMECB. This means that each 
node can receive the message from one transmitter only (and 
cannot accumulate from multiple transmitters), however one 
transmitter can transmit to multiple receivers. We had estab- 
lished in Section VII that DMECB-int is also NP complete. It 
is however interesting to note that DMECB-int needs to solve 
a weighted set cover problem when allocating powers as well; 
we know that set cover problem is o(log(n)) inapproximable 
[20], so the non-cooperative case is o(log(n)) inapproximable, 
even when ordering is provided. Greedy algorithms exist 
[20] that give 0(log(n)) approximations for the weighted 
set cover problem, and thus provide a tight polynomial time 
approximation. Therefore, to simulate a smart non-cooperative 
algorithm, we use Dijkstra's algorithm-based ordering and 
the DMECT algorithm of section VI, with the exception that 
instead of using an LP we use the greedy algorithm for power 
allocation. 

The performance comparison between our proposed cooper- 
ative algorithm and the smart non-cooperative algorithm, for 
different values of n is shown in Figure 5 and the power- 
delay tradeoff for cooperative and non-cooperative algorithms 
are presented in 6. As can be seen, the cooperative algorithm 
outperforms the non-cooperative algorithm, and the advantage 
is more pronounced when a delay constraint is imposed. 

The performance gains obtained by using MIA is shown in 
Figure 7, for a sample network of 30 nodes. 

We next study the power-delay tradeoff of the cooperative 
algorithm for different channel conditions and different values 
of network density p (in nodes/area). Figure 8 and Figure 
9, show results for EA and MIA, respectively. These figures 
highlight the sensitivity of the dense networks and those 
with poor channel conditions to delay constraints and the 
importance of having smart algorithms to minimize the energy 
consumption. 
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Fig. 5: Effect of cooperation in broadcast. 



Fig. 8: Effect of network density on power-delay tradeoff. 




Fig. 6: Power-delay tradeoff in cooperative vs non-cooperative 

Fig. 9: Different ri values for information accumulation, 
case. b 1 
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Fig. 7: Energy accumulation vs mutual information accumu- 
lation. 



X. CONCLUSIONS 

We have formulated in this work the novel problem of 
delay constrained minimum energy cooperative transmission 
(DMECT) in memoryless wireless networks, encompassing 
both EA and MIA. We have shown that this problem is 
o(log n) inapproximable in broadcast and multicast cases and 
is NP-complete in the unicast case when mutual information 
accumulation is used. For the broadcast case with EA, we 
have presented an approximation algorithm. We also devel- 
oped a polynomial-time algorithm that can solve this NP- 
hard problem optimally for a fixed transmission ordering. Our 



empirical results suggest that for practical settings, a near- 
optimal ordering can be obtained by using Dijkstra's shortest 
path algorithm. We have further showed that the unicast case 
can be solved optimally and in polynomial time when EA 
is used. We have studied the energy-delay tradeoffs and the 
performance gain of MIA using simulations, and evaluated the 
performance of our algorithm under varying conditions. 

The summary of the algorithmic results developed in this 
paper are presented in Tables I and II. Besides tackling the 
unsolved problems signified by empty slots in this table, there 
are a number of other interesting directions for future work. 
In this paper we have focused on the static problem with 
full channel-state information, which allows for centralized 
decision making. Distributed solutions to this problem which 
would be particularly suitable for more dynamic settings. 
From an analytical perspective, deriving inapproximability 
results for the NP-complete cases in the tables, as well as 
tighter approximation results are interesting avenues for further 
research. Evaluating the proposed algorithms under more re- 
alistic settings (through more detailed simulations of physical 
layer implementation or through direct implementation on 
software radio platforms, and the use of more realistic energy 
models) would certainly help in moving this work towards 
practice. Finally, our work, like most work in this domain of 
cooperative broadcasts, has focused on the single flow setting. 
It is of interest to study generalizations that allow for multiple 
simultaneous flows in the network. 
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